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Abstract 

In this paper we summarize loop quantum gravity (LQG) and we show how ideas developed in 
LQG can solve the black hole singularity problem when applied to a minisuperspace model. 

Introduction 

Quantum gravity is the theory by which we try to reconcile general relativity and quantum mechanics. 
Because in general relativity the space-time is dynamical, it is not possible to study other interactions 
on a fixed background. The theory called "loop quantum gravity" (LQG) [T is the most widespread 
nowadays and it is one of the non perturbative and background independent approaches to quantum 
gravity (another non perturbative approach to quantum gravity is called asymptotic safety quantum 
gravity [5]). LQG is a quantum geometrical fundamental theory that reconciles general relativity and 
quantum mechanics at the Planck scale. The main problem nowadays is to connect this fundamental 
theory with standard model of particle physics and in particular with the effective quantum field 
theory. In the last two years great progresses has been done to connect LQG with the low energy 
physics by the general boundary approach [3], [4]. Using this formalism it has been possible to 
calculate the graviton propagator in four [5] , [5] and three dimensions [7] . In three dimensions it has 
been showed that a noncommuative field theory can be obtained from spinfoam models [5]. Similar 
efforts in four dimension are going in progress [5]- Algebraic quantum gravity, a theory inspired by 
LQG, contains quantum field theory on curved space-time as low energy limit |10j . About an unified 
theory of particle physics and gravity authors in a recent paper have showed that spinfoam models 
[lj, including loop quantum gravity, are also unified theories, in which matter degrees of freedom are 
automatically included and a complete classification of the standard model spectrum is realized. 

Early universe and black holes are other interesting places for testing the validity of LQG. In the 
past years applications of LQG ideas to minisuperspace models lead to some interesting results in 
those fields. In particular it has been showed in cosmology [12], [13] and recently in black hole physics 
[14] , [15] , [16] , [17] that it is possible to solve the cosmological singularity problem and the black hole 
singularity problem by using tools and ideas developed in full LQG. Other recent results concern a 
semiclassical analysis of the black hole interior [TH] and the evaporation process [TU] . 

We can summarize the "loop quantum gravity program" in the following research lines 

• the first one dedicated to obtain quantum field theory from the fundamental theory rigorously 
defined; 



the second one dedicated to apply LQG to cosmology and black holes where extreme energy 
conditions need to know a quantum gravity theory. 



The paper is organized as follow. In the first section we briefly recall loop quantum gravity at 
kinematical and dynamical level. In the second section we recall a simplified model [14j showing how 
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quantum gravity solves the black hole singularity problem. In the third section we summarize "loop 
quantum black hole" (LQBH) |17j . This is a minisuperspace model inspired by LQG where we quantize 
the Kantowski-Sachs space-time without approximations. This model is useful to understand the black 
hole physics near the singularity because the space-time inside the event horizon is of Kantowski-Sachs 
type. 



1 Loop quantum gravity in a nutshell 

In this section we recall the structure of the theory introducing the Ashtekar's formulation of general 
relativity [20] , the kinematical Hilbert space, quantum geometry and quantum dynamics. 



1.1 Canonical gravity in Ashtekar variables 

The classical starting point of LQG [T] is the Hamiltonian formulation of general relativity. In ADM 
Hamiltonian formulation of the Einstein theory the fundamental variables are the three-metric q a b of 
the spatial section £ of a foliation of the four-dimensional manifold A4 = M x E, and the extrinsic 
curvature K a b- In LQG the fundamental variables are the Ashtekar variables: they consist on an 
SU{2) connection A l a and the electric field Ef, where a, b, c, ■ •■ — 1, 2, 3 are tensorial indices on the 
spatial section and i, j, k, ■ ■ ■ = 1, 2, 3 are indices in the su(2) algebra. The density weighted triad Ef 
is related to the triad e l a by the relation Ef = ^e abc eijk e\ e k c . The metric is related to the triad by 
q a b = e a e t s ij ■ Equivalently, 

y/d^q)q ab = E° E)^. (1) 
The rest of the relation between the variables (A l a ,Ef) and the ADM variables (q a b,K a b) is given by 

Ai = Ti + 7 K ab E b j 5^ (2) 



where 7 is the Immirzi parameter and T l a is the spin connection of the triad, namely the solution of 

a i r A yj k 
-f>] + 6 jk 1 [a C b] 



Cartan's equation: <9[ a e?, + e) k P? e?, = 0. 



The action is 



K7 



J dt J d 3 x -2Tr(£ a i a ) -NH- N a H a - N*Gi , (3) 



where N a is the shift vector, N is the lapse function and N l is the Lagrange multiplier for the Gauss 
constraint Gi- We have introduced also the notation E = E a d a = Efr l d a and A — A a dx a = A l a T % dx a . 
The functions Ti, T~L a and Qi are respectively the Hamiltonian, diffeomorphism and Gauss constraints, 
given by 

H{Ef,A\) = -4 e- 1 Tr {F ab E a E b ) - 2 e" 1 (1 + 7 2 ) EfE b K\ a K\ 



b] 

7~lb{Ef, A* ) = Ej F° ab - (1 + j 2 )K i a G i 

g i (E?,Ai) = d a E? + e>? j AiE%, (4) 

where the curvature field strength is F a b = d a Ab — dbA a + [A a , A ] and e = det(e^). The constraints (JH) 
are respectively generators for the foliation reparametrization, for the S surface reparametrization and 
for the gauge transformations. The symplectic structure for the Ashtekar Hamiltonian formulation of 
general relativity is 

{E^(x),Al(y)} = K-ySSfyfay), {E*(x), E b (y)} = {A{(x), A\{y)} = 0. (5) 

General relativity in metric formulation is defined by the Einstein equations = 8itGn T^ u . The 
Ashtekar Hamiltonian formulation of general relativity is instead defined by the constraints H = 0, 
H a = 0, Gi = and by the Hamilton equations of motion: A a — {A^,7i} and Ef = {Ef,H}. 
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1.2 The Dirac program for quantum gravity 



The general strategy to quantize a system with constraints was introduced by Dirac. The program 
consist on : 

1. find a representation of the phase space variables of the theory as operators in an auxiliary 
kincmatical Hilbert space Hkin satisfying the standard commutation relations, i.e., { , } — > 
-i/h[ , ]; 

2. promote the constraints to (self-adjoint) operators in Hkin- For gravity we must quantize a set 
of seven constraints Qi(A, E), H a {A, E), and H(A, E) and we must solve the quantum Einstein's 
equations (for 7 = i) 

H\il>) = [-Ae- 1 TT(F ab E a E b )+n M ] W = 0, 



Mb 



W) = 0, 



dW)= [d a Ef + 4 AiE a k + g Mi ] |V) =0 (6) 
We will consider pure gravity then the matter constraints are identically zero. 
3. introduce an inner product defining the physical Hilbert space H p h ys - 



1.3 Kinematical Hilbert space 

The fundamental ingredient of LQG is the holonomy of the Ashtekar connection along a path e, 
h e [A] = Pexp — J e A E SU(2). Given two oriented paths e\ and such that the end point of e\ 
coincides with the starting point of so that we can define e = we have the composition rule 
h e [A] = h ei [A]h e2 [A]. By a gauge transformation the holonomy transforms as 

h' e [A]=g{x{Q)) h e [A] 9-\x{\)), (7) 

and by a Diffeomorfism of the three dimensional variety <fi G Diff(S) we have 

h e [<f?A] = ht-i [e) [A], (8) 

where <f>*A denotes the action of <f> on the connection. In other words, transforming the connection 
with a diffcomorphism is equivalent to simply moving the path with <j> . 

We introduce now the space of cylindrical functions (Cyl 7 ) where 7 denotes a general graph. 
A graph 7 is defined as a collection of paths e C S (e stands for edge) meeting at most at their 
cndpoints. If N e is the number of paths or edges of the graph and ej, for i = 1, • • • N e , are the edges 
of the corresponding graph 7 a cylindrical function is an application / : SU(2) Ne — > C, defined by 

V 7l/ L4] := f(h ei [A], h e2 [A], - ■ ■ h eNc [A]). (9) 

Two particular examples of cylindrical functions are the holonomy around a loop, VK-J^] := Tr[/i 7 L4]] 
and the three edges function Ql'^e, i A \ - \A]f ] 1/2 @(h e2 [A]) AB ^ 2 &(h e3 [A]) CD f£ BCD , 

where ^(h ei ) is the SU (2) representation for the holonomy along the path ei and f ABCD are complex 
coefficients. The algebra of generalized connections is given by Cyl = U 7 Cyl 7 . 

We introduce now the space of spin networks states. We label the set of edges e C 7 with spins 
{j e }- To each node n C 7 one assigns an invariant tensor, called intertwiner, i n , in the tensor product 
of representations labelling the edges converging at the corresponding node. The spin network function 
is defined 

*r,{i.},K}[4 =0^(g) j °nhe[A\) , (10) 

nC7 eC7 
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where the indices of representation matrices and invariant tensors are implicit to simplify the notation. 
An example of spin network state is 

[A] = l @(h ei ^ 2 ^(h e2 [A]) AB l '^{K 3 [A]) C D °t C <rf D , (11) 

where for the particular representations converging to the two three-valent nodes of the graph the 
intertwiner tensor is the Pauli matrix. The spin network states are gauge invariant because for any 
node of the graph we have invariant tensors (intertwiners), then on the spin network states the Gauss 
constraint is solved as asked from the Dirac program of the previous subsection. 

To complete the Hilbert space definition we must introduce an inner product. The scalar product 
is defined by the Ashtekar-Lewandowski measure 

<ip~t,f,ip~t>,g >= MA£C0-y,/Vy, ff ) =/ II dh e f( h e 1 ,---h eN Jg(h ei ,...h eN J, (12) 

err , 

where we use Dirac notation and f(h ei ,...h eN ), g(h ei , . . . h eN ) are cylindrical functions; T 77 / is 
any graph such that both 7 C T 77 / and 7' C T 77 /; dh e is the Haar measure of SU(2). The scalar 
product (|12[) is non zero only if the two cylindrical functions have support on the same graph. The 
kincmatical Hilbert space Hkin is the Cauchy completion of the space of cylindrical functions Cyl in 
the Ashtekar-Lewandowski measure. In other words, in addition to cylindrical functions we add to 
Hkin the limits of all the Cauchy convergent sequences in the norm defined by the inner product : 

i' = En=i*mPn, IMI 2 = £^iKl a IWI 2 - 

We complete the construction of the theory at kinematical level solving the diffcomorphism con- 
straint. Given ip-yj € Cyl the finite action of a Diff. transformation is implemented by an unitary 
operator Ux> such that 

Wj,[# 7l/ [i]=^- v [A]. (13) 

The states invariant under Diff. transformations satisfy Uj}[(j)]ijj = 4> and are distributional states in 
the dual space of Hunt "0 G Cyl* 

([V> 7 ,/]l= E <^fPv[4>]= £ <Vwl, (14) 

06Diff(S) 0SDiff(S) 

where the sum is over all diffeomorphisms which modified the spin network. The brackets in | de- 

note that the distributional state depends only on the equivalence class [tp-yj] under diffeomorphisms. 
Clearly we have ([^ 7 j]\U v [a] = ([ip 7 ,f]\ V«£ Diff(E). 

We conclude that the Dirac's program at kinematical level is realized by the Gelfand triple 

Cyl C H kin C Cyl* 5 ^ 2) Cyl C H km C Cyl^ D -^' H ms C Cyl*, (15) 

where Cylo is the subspace of cylindrical functions invariant under SU(2). 

At quantum level the phase space variables operators are represented on the spin network space 
by the holonomy operator h e [A] that acts multiplicatively on the states and by the smearing of the 
triad Ef on a two dimensional surface S £ E 

E[S,a] = J da^^^E t e ahc = -itery J ^da^^a^e abc , (16) 
5 s 

where a 1 is the smearing function with values on the Lie algebra of SU(2). The action of E[S, a] on 
the spin network states can be calculated using 



E[S,a]h e (A) — —il%'-fa l h ei (A)Tih e2 (A), e = e\ U 
5 

'6Ai 



E?(x)h e (A) = - t l 2 p J-h e (A) = I d S e a ( S )<5 3 (e( S ),x)/ lei (A)r^ e2 (A), (17) 
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Figure 1: In this picture we show the spin network physical meaning. The graph on the left represents a 
particular spin network. In the center we represent the same spin network and the dual decomposition 
of the space section in chunk of space. In the last picture on the right we consider another spin network 
and a particular dual volume region. The yellow region is a chunk of space with volume eigenvalues 
related to the red intertwiners and area eigenvalues given by the SU(2) representations associated to 
the green edges. 

where e a (s) is tangent to the curve e(s) in the graph 7. The pair (/i e [A], E[S, a]) realizes the first 
point of the Dirac's program. 

1.4 Quantum geometry and dynamics 

We are going to give a physical interpretation of the Hilbert space previously introduced. We consider 
the spatial section S of the space-time and we study the spectrum of the area S and volume R in the 
section £ 21 . We define the area of a surface S as the limit of a Riemann sum 

N 

A s = lim Ag , A% = J2 VMSiWW) (18) 

7=1 

where N is the number of cells. The quantum area operator is As = limjv^oo ^.5^ Using (|17|) we 
calculate the action E l (S I )E i (Si) j &{h e [A]) mn = (l 2 P j) 2 (j{j + 1)) j @{h e [A]) mn . The area spectrum 
for spin network without edges and nodes on the surface is 

AsllJe^n) = l 2 pl^2\Jjptip + l)l7>Je,tn), (19) 

where j p are the representations on the edges that cross the surface S. Now we consider a region R 
with a number n of nodes inside. The spectrum of the volume operator for the region R is 

VR|7,ie,0 = lW 2 E \J™ M Un)hJe,ln). (20) 
nCR 

We have all the ingredients to give a physical interpretation of the Hilbert space. The states have 
support on graphs that are a collection of nodes and edges converging in the nodes. The dual of a 
spin network corresponds to a discretization of the three dimensional surface S. The dual of a set of 
edges is the 2-dimensional surface crossed by the links and the dual of a set of nodes is the volume 
chunk contained nodes (see Fig[T]). 

We must now implement the Hamiltonian constraint on the Hilbert space. The Euclidean part of 
the constraint S E (N) = J^d 3 xN(x)H(Ef, A 4 Q ), 

s!is, f'" w H'" A (2i » 
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Using the Thiemann's trick pQ we can express the inverse of y det(E) by 



E b E c A 



and the Hamiltonian constraint is 



^ (iV ) = — dx 3 N(x) t ahc ^Fl h {Ai, V) . 
K l J 



(22) 



(23) 



Now we define the Hamiltonian constraint in terms of holonomies. Given an infinitesimal loop a a b on 
the a6-plane in the surface £ with coordinate area e 2 , we can define F* b in terms of holonomies by 
K ab [A] - h~l b [A] = e 2 F* ab T t + C(e 4 ) and h£[A] {h Ca [A], V} = e {A a , V) + 0{e 2 ) (e a is a path along 
the a-coordinate of coordinate length e) . With these ingredients the quantum constraint can formally 
be written 



S E (N) = — lim V Nj e afcc Tr \(h aI [A]~hJ [A]jhJ[A] \h e i[A],V 

e— >0 t-r^ L ab ab c L 



(24) 



where we have replaced the integral by a Ricmann sum over cells of coordinate volume e 3 . It is easy 
to see that the regularized quantum scalar constraint acts only on spin network nodes, because in (|2"T]) 
-F a b and E a E b are respectively antisymmetric and symmetric in indexes on spin network states. In 
fact E a E br ^ 1 j ~ e a e b 4>j'j' (this is a consequence of (|T7)l ). The action of (|2"4"|) on spin network states 
is S e (N)i/)yj — J2 n ~/N n S™ tp-yj, where S 1 ™ acts only on the node n C 7 and N n is the value of the 
lapse N(x) at the node. The scalar constraint modifies spin networks by creating new exceptional 
links around nodes. The Euclidean constraint action on 4-valent nodes is [T] 



ST 



jklm 7 opq 




/ J Sjmkl,opq 



(25) 



The only dependence on the regularization parameter e is in the position of the extra edges in the 
resulting spin network states, then the limit e — > can be defined on diffcomorphism invariant states 
in Hjjiff. The key property is that in the diffcomorphism invariant context the position of the new 
link is irrelevant. Therefore, given a diffeomorphism invariant state ([(f>]\ € H^iff C Cyl*, the quantity 
(<j>\S e (N)\ip) is well defined and independent of e. 

The operator (f2~4"|) defines the dynamics and an unitary implementation of the constraint S E (N) 
gives the evolution amplitude from a spin network s to a new spin network s'. Introducing the 
projector operator P = J D[N) exp (i J N(x)S E (x)) we can characterizes the solutions of quantum 
Einstein equations by P\s), V|s) € Pikm- The matrix elements of P define the physical scalar product, 
W(s,s') := (s|s') phys '■— { S \P\ S ') ■ The amplitude W(s,s') solves the Hamiltonian constraint in the 
following sense. If \^ P h ys ) ■= P\s') we have that S E (N)\^ phys ) = 0, but {s\9 phya ) := 9 p hy a (s) = 
W(s, s'), then we obtain 



(s\S E (N)\* phya ) =^( s |5 £ (iV)| s ")( 5 "l* P ^> = 

s" 

= S E (N) ss „y phys ( s ") = S E (N) ss „W(s", s') 



(26) 



Relation ([2"6"|) shows the amplitude W{s, s') is in the Hamiltonian constraint kernel. W{s, s') realizes 
the Dirac's program because corresponds to the finite implementation of the scalar constraint on the 
kinematical states and defines a class of models called spinfoam models [1]. 
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2 Avoidance black hole singularity in quantum gravity 



In this section we study the black hole system inside the event horizon in ADM variables considering 
a simplified minisuperspace model [2] and in particular using the fundamental ideas suggested by 
full LQG and introduced in the first section. The simplification consist on a semiclassical condition 
which reduce the phase space from four to two dimensions. This is an approximate model but it is 
useful to understand the ideas before to quantize the Kantowski-Sachs system in Ashtekar varables. 

2.1 Classical theory 

Consider the Schwarzschild solution inside the event horizon; the metric is homogeneous and it reads 

dg2 = - (M/jT i) + i^r 1 l ) dr ' 2 + T2 ( sin2 M ^ 2 + d ° 2 )- ( 2? ) 

This metric is a particular representative of the Kantowski-Sachs class [25 

ds 2 = -clt 2 + a{t) 2 dr 2 + &0) 2 (sin 2 9d4> 2 + d6 2 ). (28) 

Introducing (|28p in the Einstein-Hilber action we obtain S = — 2 q J dt ab 2 + 2dbb — a (where 

R = J dr is a cut-off on the radial cell) [H] . We introduce in the action the classical relation 
a 2 = 2MG N /b(t) - 1, obtaining 



R 
2G~ 



N 



dt 



Vb ( x b yh - 2 ymu^ 



y/2MG N V 2MG N ) ^b V 2MG N 

-1/2 



(29) 



The momentum conjugates to the variable b(t) is p = G ^/2MG — 2MG N ) ^' an< ^ ^ ne Hamil- 
tonian constraint, by Legendre transform, is 



H =pb- L 



G NP 2 R 



2R 2G 



N 



Vb V 2MG 



N 



1/2' 



(30) 



We introduce a further approximation. In quantum theory, we will be interest in the region of the scale 
the Planck length lp around the singularity. We assume that the Schwarzschild radius r s — 2MGn 
is much larger than this scale. In this approximation we can write 1 — b/2MGN ~ 1 and H becomes 



H = 



G N p 2 R \ V2MG 



2R 2Gn J \/b 



In the same approximation the volume of the space section is V — 4:TtR\/2MGn b 3 ^ 2 :— l 6 3 / 2 . 

The canonical pair is given by b = x and p, with Poisson brackets {x,p} = 1. We now assume 
that x G M. This choice it not correct classically, because for x — we have the singularity, but it 
allows us to open the possibility that the situation be different in the quantum theory. We introduce 
an algebra of classical observables and we write the quantities of physical interest in terms of those 
variables. We are motivated by loop quantum gravity to use the fundamental variables x and 

U s (p) = cxp( — f-ip) (31) 

where 6 is a real parameter (see next paragraph for a rigorous mathematical definition of 6) and L 
fixes the unit of length. The operator ([3"T|) can be seen as the analog of the holonomy variable of loop 
quantum gravity. 
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A straightforward calculation gives 



{x,U s (p)} = 8TrG N l -^U s (p), 

U^{V n , U S } = q U^{\x\^,U s } = i 8t:G n Z " —aga^x]*- 1 . (32) 

These formulas allow us to express inverse powers of x in terms of a Poisson bracket between Us and 
the volume V, following Thiemann's trick [26] ■ For n = 1/3 ([3"2"| gives 

1M TmW"' { ' %) - (33) 

We will use this relation in quantum mechanics to define the physical operators. We are interested to 
the quantity l/|x| because classically this quantity diverge for |x| — > and produce the singularity. 
We are also interested to the Hamiltonian constraint and the dynamics and we will use l|33p for writing 
the Hamiltonian. 



2.2 Polymer quantization 

In this section we recall the polymer representation [22] of the Weyl-Heisenberg algebra and we com- 
pare this representation with full LQG. The polymer representation of the Weyl-Heisenberg algebra is 
unitarily inequivalent to the Schroedi-nger representation. Now we construct the Hilbert space Hpoly 
First of all we define a graph 7 as a finite number of points {xi} on the real line R. We denote by 
Cyl 7 the vector space of function f(k) (/ : R — ► C) of the type 

/(*) ,rk (34) 

i 

where k € R, Xj € R and fj £ C and j runs over a finite number of integer (labelling the points of 
the graph). We will call cylindrical with respect to the graph 7 the function f(k) in Cyl 7 . The real 
number k is the analog of the connections in loop quantum gravity and the plane wave e~ tkx i can be 
thought as the holonomy of the connection k along the graph {xj}. 

Now we consider all the possible graphs (the points and their number can vary from a graph to 
another) and we denote Cyl the infinite dimensional vector space of functions cylindrical with respect 
to some graph: Cyl = M GyLy. A basis in Cyl is given by e~ lkXj with ^ e -" lkx i\ e ~' lkx j''j = S XiiXj . Hp a i y 
is the Cauchy completion of Cyl or more succinctly Hp a i y = L2(M.BohridfJ,o), where Rs ; lr is the 
Bohr-compactification of R and dfio is the Haar measure on Rsohr- 

The Weyl-Heisenberg algebra is represented on Hp Q iy by the two unitary operators 

V(X)f(k) = f(k-X), U(S)f(k)=e M f(k), (35) 

where A, S € R. In terms of eigenkets of V(X) (we associate a ket \xj) with the basis elements e~ lkXj ) 
we obtain 

V(X)\ Xj ) = e iX *i\ Xj ) , U{5)\x 3 ) = \ Xj - 6). (36) 

It is easy to verify that V(X) is weakly continuous in A, whence exists a self-adjoint operator x such 
that x\xj) = Xj\xj) gg, [21]. 

The operator analogy between loop quantum gravity and polymer representation is the following: 
the basic operator of loop quantum gravity, holonomies and electric field fluxes, are respectively 
analogous to the operators U(5) and x with commutator [x,U(5)} = —5U(S). The commutator is 
parallel to the commutator between electric fields and holonomies. As, in the polymer representation, 
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the unitary operator U(S) is well denned but the operator p doesn't exist, in loop quantum gravity 
the holonomies operators are unitary represented self-adjoint operators but the connection operator 
doesn't exist. As x, the electric flux operators are unbounded self-adjoint operators with discrete 
eigenvalues. 

After this short review on polymer representation of the Weyl-Heisenberg algebra we return to our 
system. 



2.3 Polymer black hole quantization 

Following the previous section we quantize the Hamiltonian constraint and the inverse volume operator 
in the Polymer representation of the Weyl-Heisenberg algebra. The operators are x, acting on the 
basis states according to 

x\(j) = Lfi\fi) , (n\v) = (37) 

(we have redefined the continuum eigenvalues of the position operator of the previous section xi — > fi) 
and the operator corresponding to the classical momentum function Us — e l87rGjv<5p / i . We define the 
action of Us on the basis states using the definition (|37[) and using a quantum analog of the Poisson 
bracket between x and Us 

Us\fx) = \fx-S) , [x, Us] = -SLUg. (38) 

Using the standard quantization procedure [, ] — > ih{ , }, the Poisson bracket (|32|) and (l38|) we obtain 
the value of the length scale L = ^/%tFGnH = lp. 



2.3.1 Avoidance black hole singularity and regular dynamics 

We recall that the dynamics is all in the function b(t), which is equal to the the radial Schwarzschild 
coordinate inside the horizon. The important point is that b(t = 0) = and this is the Schwarzschild 
singularity. We now show that the spectrum of the operator does not diverge in quantum me- 
chanics and therefore there is no singularity in the quantum theory. 

Using the relation (l3"3"l) . and promoting the Poisson brackets to commutators, we obtain (for 8—1) 
the operator 



1 1 / • r , 

M 2irG N hll 



(u~ l V^,U Y ■ (39) 



The action of this operator on the basis states is (the volume operator is diagonal on the basis states, 
V-| / i)=ioWt| / i)=io|£rf«|rt) 



HW-te^-i'-'i'l'w- (40) 

We can now see that the spectrum is bounded from below and so we have not singularity in the 
quantum theory. In fact the curvature invariant TZ flvpcr TZ tl ' yp ' 7 = 48M 2 G%/x(t) 6 is finite in quantum 
mechanics in \i = 0. The eigenvalue of the operator \j\x\ for the state |0) corresponds to the classical 
singularity and in the quantum case it is 4/l p , which is the largest possible eigenvalue. For this 
particular value the curvature invariant it is not infinity 

K^WIO) = |a;|6 " 10) = ^ N \0). (41) 
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If we consider the h — > limit we obtain the classical singularity so the result is a genuinely quantum 
gravity effect). On the other hand, for \fi\ —> oo the eigenvalues go to zero, which is the expected 
behavior of l/\x\ for large \x\. 

Now we study the quantization of the Hamiltonian constraint near the singularity, in the approx- 
imation (|2.ip . There is no operator p in polymer quantum representation that we have chosen, hence 
we choose the following alternative representation for p 2 . Consider the classical expression 



-, 7z — 777 lim 

{8irG N ) 2 s^o 



2 -Us- Ug 
S 2 



(42) 



We have can give a physical interpretation to 8 as 8 = l v j 'L p hy S , where Lph ys is the characteristic size 
of the system. Using (|4"2"j) we write the Hamiltonian constraint as 



H 



C 



^ /3 Li/2 



Us + tig 1 - (2-C)!] sgn(a;) (if- 1 



(43) 



where C 



■5/2 , 



and C 



The action of H on the basis states is 



H\») = C V( M ) [\h -S) + |/i + 8) - (2 - C')|M>] , 

for u ^ 



V( M ) 



|| M _ 5 |1/2_| M |1/2| 
1/2 



for pi = 



(44) 



We now calculate the solutions of the the Hamiltonian constraint. The solutions are in the dual 
space of Hp iy. A generic element of this space is = ip(p)(p\. The constraint equation 

H\ip) = is now interpreted as an equation in the dual space {ip\W ; from this equation we can derive 
a relation for the coefficients V'O-O 



V(/i + 8) ip(fi + S) + V(ji - 8) V(M -8) -(2- C) V(p) ipifJ.) = 0. 



(45) 



This relation determines the coefficients for the physical dual state. We can interpret this states as 
describing the quantum spacetime near the singularity. From the difference equation (|45|1 we obtain 
physical states as combinations of a countable number of components of the form ip(fi + nS)\fi + n8) 
(8 ~ Ip/Lphys ~ 1); any component corresponds to a particular value of volume, so we can interpret 
ip(/j, + 8) as the wave function describing the black hole near the singularity at the time u + 8. A 
solution of the Hamiltonian constraint corresponds to a linear combination of black hole states for 
particular values of the volume or equivalently particular values of the time. 



3 Loop quantum black hole 

In this section we quantize the Kantowski-Sachs space-time in Ashtekar variables and without approx- 
imations [T7] . 

3.1 Ashtekar variables for Kantowski-Sachs space-time 

The Kantowski-Sachs space-time is a simplified version of an homogeneous but anisotropic spacetime, 
written in coordinates (t,r,9,<f>). An homogeneous but anisotropic space-time of spatial section S of 
topology S = K x S 2 is characterized by an invariant connection 1-form of the form [57], [IE] 

A {1] = A r (t)T S dr+(A 1 (t)n+ A 2 (t)T 2 )de+(A 1 (t)T 2 - A 2 (t)n) sin9 # + r 3 cos 6 d<j). (46) 
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The Ti are the generators of the SU(2) fundamental representation. They are related to the Pauli <Ji 
matrix by Ti = —\<Ji- On the other side the dual invariant densitized triad is 



E m =E r (t)T 3 sine— + {E\t)r l +E 2 {t)T 2 ) sin0 — + {E l {t) r 2 - E 2 (t)n)— . (47) 
or ou 



Since spacetime is homogeneous, the diffeomorphism constraint is automatically satisfied. 

In this paper we study the Kantowski-Sachs space-time with space section of topology R x S , 
the connection Am is more simple than in (|4"6")l with A 2 = Ax, and in the triad (|4T|) we can choose 
the gauge E 2 = E 1 29J . There is a residual gauge freedom on the pair (AxjE 1 ). This is a discrete 
transformation P : (Ai : E r ) — > (— Ax,— E 1 ); we have to fix this symmetry on the Hilbert space. 
The Gauss constraint is automatically satisfied and the Euclidean part of the Hamiltonian constraint 
becomes 



2\/2sin6» sgn(£ r 

He — 



2A r E r A 1 E 1 + (2(A 1 ) 2 - l)^ 1 ) 2 (48) 



wc redefine E r = E and A r = A. The connection between the metric (|28|) and the density triad is 
q a b = d\&g{2{E 1 ) 2 /\E\, \E\, \E\ sin 2 9). Another useful quantity is the volume of the spatial section £ 

V = J drd<pd9 ^/det(q) = 4:TrV2R^/\E\\E 1 \, (49) 

where R is a cut-off on the space radial coordinate. The spatial homogeneity enable us to fix a linear 
radial cell C r and restrict all integrations to this cell [12] • To simplify notations we restrict the linear 
radial cell to the Planck length lp and so we can take J dr — R = lp in the action functional (|4"9l . 
Now we are going to use R= lp in all the paper. 

The classical symplectic structure of the phase space can be obtained by inserting the symmetry 
reduced variables in the action This gives 

s = l- 

«7 



— J dt J drd(/)dO sind Tr [-2E a A^ + . . . = — ^ J dt EA + 4E 1 A 1 + . . . . (50) 



We can read the symplectic structure of the classical phase space directly from the reduced action 
(|50|) . The phase space consists of two canonical pairs A,E and A\,E and from ([5U|) we can obtain 
the simplectic structure. The simplectic structure is given by the poisson brackets, {A, E} = 
and {AijE 1 } = . The coordinates and the momenta have dimensions: [A] = [Ai] = L°, 

[E] = L 2 and [E 1 ] =L. 

The elementary configuration variables used in LQG are given by the holonomics along curves in 
the spatial section Kx5 2 and the fluxes of triads on a two-surface inRxfi* 2 . We restrict our attention 
to three sets of curves. More precisely, we consider only spin networks [T] based on graph made just 
of radial edges, and of edges along circles in the ^-direction or at 9 — tt/2. 

Let us introduce the fiducial triad °ef — diag(l, 1, sin -1 9) and co-triad °uj I a — diag(l, 1, sin#). The 
holonomy along a curve in the direction "J" is given by hi = exp J A l jTi, 

hi = expLA/^p r 3 ] , h 2 = exp[Ai^ fa + n)] , h 3 = exp[Ai/x (r 2 - n)] , (51) 

where A[ = (0, 0, A), A\ = (Ai,Ai,0) and A\ = {—Ax, Ax, 0). The connection in ((SU is integrated in 
the direction "7"; fiolp is the length of the curve along the direction r, is the length of the curve 
along the directions 9 and <f>. The length are defined using the fiducial triad °ef . 

Recall that the Hamiltonian constraint can be written in terms of the curvature F a j, and the Poisson 
bracket between A a and the volume V 26 1 . The Euclidean part of the Hamiltonian constraint becomes 

H E = -— [ d 3 xNe abc Tr[F ab {A c ,V}]. (52) 
K 1 J 
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Because of homogeneity we can assume that the lapse function N is constant and in the rest of the 
paper we will set N = 1. 

We can express the curvature F a b and the Poisson bracket r\ {A^, V} in terms of holonomies [T7] 
obtaining the following form for the Hamiltonian constraint 

H E = eIJKTr [hihjhfhfh^j] h^{h K ,V}] , (53) 

K7 ^o UK 

where hyj] = exp(-/i^ C/jt 3 ) and C u = S 2 iS 3J - 6 3l 6 2 j- 
Using the classical identity 

JXTT^v = ~ vT: lp) l <*> ^ viH4, vhiA*, viy, (54) 

(det(g a b)j4 3! K d 7 d (smy)2 
we can define the inverse of the volume in terms of holonomies. The result is |17j 



1TT7^ = T?Tp% E ^ K ^hfih;, V*}] T^h-/ {hj , Vi}] Tr[rV{h K , V*}]. 
(det(g))4 3K 3 7 3 ^(sm6»)2 fj£ 

(55) 

The definition (l55|) will be useful to calculate the inverse volume spectrum. 



3.2 Quantum Theory 

We construct the kinematical Hilbert space H^m for the Kantowski-Sachs minisuperspace model in 
analogy with the full theory. As in the second section we define a graph T as a finite number of couple 
of points (fiEi j A*© 1 *)) where fiEi, ^E 1 i € ®- We denote by Cyl r the vector space of function f(A, A\) 
(/ : R 2 -► C) of the type 

/(A A) = £ fa e~^ + ^^. (56) 

where A,A\ € R, , ^E 1 j G R i /y € C and i,j run over a finite number of integers (labelling 
the points of the graph). We call the function f(A,A\) in Cylr cylindrical with respect to the graph 
r. We consider all possible graphs (the points and their number can vary from a graph to another) 
and denote by Cyl the infinite dimensional vector space of functions cylindrical with respect to some 
graph: Cyl = Ur^ylr- Thus, any element f(A,A\) of Cyl can be expanded as in (|56"|) . where the 

uncountable basis e 2 (£> e ^ is now labeled by arbitrary real numbers (he , He 1 )- A basis 
in Cyl is given by IhejI^e 1 ) = \he) <8> \he 1 )- Introducing the standard bra-ket notation we can define 
a basis [13] in the Hilbert space via 

{A\he) <8> {Ai\fi B i) = e ^ E 2 pA ^e-^, (57) 

The basis states (|57|) are normalizable in contrast to the standard quantum mechanical representation 
and they satisfy 

(fiE,fJ>E 1 \ l/ B,VE 1 ) =8 flB ,u B 5 flEl ,u Bl - (58) 

The Hilbert space Tikin is the Cauchy completion of Cyl or more succinctly Hkin = L2(R 2 Bohr , d/io), 
where M-Bohr is the Bohr-compactification of K and d^o is the Haar measure on R% ohr - In LQG (or 
in polymer representation) the fundamental operators are E, E , hj and 
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• the momentum operators can be represented on the Hilbcrt space by 

A -lip d pi .lip d 

E ~* 77 i E — * ~ l ~T7. TJ-i (59) 

4tt dA 16tt dAi 

and the he spectrum of these two momentum operators on the Hilbert space basis is 

E\fi E ,(iE 1 ) = — 5 Imb,Mbi), E \ii E ,fi E i) = — — =- (60) 

07T 167TV2 



• the holonomy operators in the directions r, 0, of the space section lx5 2 are : h^ E \ h^ 1 

and h^ 1 , where fJ. E lp is the length along the radial direction r and fi El is the length along 
the directions 8 and </> (all the length are define using the fiducial triad °e"). The holonomies 
operators act on the Hilbert space Jikin by multiplication. 

We have to fix the residual gauge freedom on the Hilbert space. We consider the operator P : 
\HEi A'b 1 ) ~~ > |MBj — ^e 1 ) an d we impose that only the invariant states (under P) are in the kinematical 
Hilbert space. The states in the Hilbert space are : [|mb, Mb 1 ) + 1Mb, "Mb 1 )] f° r Mb 1 7^ an d the 
states |/i£,0) for fx E i = 0. 



3.2.1 Inverse volume spectrum 



In this section we study the black hole singularity problem in loop quantum gravity calculating the 
spectrum of the inverse volume operator. The operator version of the quantity sgn(E)/ (det(q))* 
defined in the formula (l55l) is 



sgn(E) 256i(4irl P )ie 



(det(?))i 3l pl 3 fi 3 {sm6)i 



UK 



UK 



Tr 



1/2" 



Tr 



r>hj l [hj,Vv] Tr T*V[fcjf,V] 



(61) 



To calculate the action of (|6"Tj) on the Hilbert space basis we Introduce the normalized vectors n\ = 
(0, 0, 1), n\ = ^(1, 1) 0), rig = lj 1, 0). Using the vectors we can rewrite the holonomies hj of 
(|5Tj) as for = exp(A/ r») = cos(Ar/2) + 2n\ t% sin(A//2), where Ai—\ = Alp[i and j4/=2 = ^4/=3 = 
AifiQi/2. The action of the multiplicative operator hi on the bases states is 

fix I — irijTi I + m}Ti 

hi\Hi) = g lMi + Mo)H 5 iMi-Mo)- (62) 

The spectrum of the inverse volume operator is 

— — ^— t \hb,Hbp) = (^zY- — TT—77I ImbI* Imb 1 !* 
(det(g))4 V7/ Zp^g(sin0)2 

Ime + Mol^ - I Mb - Mo|* 

The spectrum of the inverse volume operator is bounded above, near the classical singularity which 
is in E = or fi E = 0, and reproduces the correct classical spectrum of l/(det(g))* for large volume 
eigenvalues (Fig[2|). 



(\n E i + fio\2 - \fi E i -Mol^J Ime^e 1 )- (63) 
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Figure 2: Plot of the inverse volume spectrum for fi^ 1 ^ and V/x^- 



3.2.2 Quantum dynamics 

In this section we study the dynamics of the model solving the Hamiltonian constraint in (the dual 
space of) the Hilbert space. 

The quantum version of the Hamiltonian constraint defined in (|53[) can be obtained promoting the 
classical holonomies to operators and the poisson bracket to the commutator 



hihjhJ l hj 1 h [IJ] h] ( 1 [hK,V] . (64) 



Using the relations in (f6"2"| we can calculate the action of the Hamiltonian constraint on the Hilbert 
space basis |/i£,/x B i) and solve the Hamiltonian constraint to obtain the physical states. As in 
non-trivially constrained systems, we expect that the physical states are not normalizable in the 
kincmatical Hilbert space. However, as in the full loop quantum gravity, we again have the Gelfand 
triple Cyl C Hkin C Cyl* and the physical states will be in Cyl*, which is the algebraic dual of Cyl. 
An element of this space is (-01 = M x VVe^ei (mbiMe 1 !- The constraint equation H\ip) = is 
now interpreted as an equation in the dual space (^W; from this equation we can derive a relation 
for the coefficients ^Ve^e* 

~ a HB — 2M0:M E i — 2 Mo I^P-e — 2mo,A i e i — 2/^0 a fi.E+2fi.a,l-i E i—'2na ^•Ve+2mo,^ £; i ~ 
+ q; Me-2^o,M e i+2aio ^/Ve-2^o,M e i+2aio — a ^ E +2^ ,A l E l +2^o 7 /Ve+ 2 Mo ,Me1 + 2 Mo 

, sin(/zg/2) - cos(/i§/2) 



2 

.2 



\0fJ,E,P E l— ^Ve,£Ie 1 — 4 Mo ^P^E^e 1 ^VE^e 1 iA'e 1 + 4 «) ^Ve .Me 1 + 4 Mo ) 



-2sin(/i /2) (/3 a1e , Me1 -2 Mo ll>HB,li B i-2na + P^M,H E i+2^ ^me,M e i+2m) — 0; ( 65 ) 

where the functions a{pE, He 1 ) an d P(^e, Mb 1 ) are define by 

V,m e i : = ImbI* (1Mb 1 +/- 1 o| - 1Mb 1 -Mo|) , /?mb^ b i : = Im^KImb +Mo|^ - 1Mb — Mo|*) (66) 

How can be seen from equation (|65|) the quantization program produce a difference equation and 
imposing a boundary condition we can obtain the wave function VW-m e i f° r the black hole. We can 
interpret ^ Si)1 , as the wave function of the anisotropy "E 1 " at the time "25". It is evident from 
(|66[) that the dynamics is regular in /xe = where the classical singularity is localized. As in loop 
quantum cosmology also in this case the state 0o,o decouples from the dynamics and the quantum 
evolution does not stop at the classical singularity. The "other side" of the singularity corresponds to 
a new domain where the triad reverses its orientation. 



14 



Conclusions 



In this paper we have summarized loop quantum gravity theory and we have applied the ideas to 
study the space-time region inside the Schwarzschild black hole horizon. Because the space-time 
region inside the horizon is spatially homogeneous of Kantowski-Sachs type |25j , we have studied this 
minisuperspace model. This is an homogeneous but anisotropic minisuperspace model with spatial 
topology lx5 2 . We have analyzed the model firstly in ADM variables with some drastic simplification 
and then in Ashtekar variables. We have quantized the reduced model using a quantization procedure 
induced by the full "loop quantum gravity" . Our analysis it was useful in order to understand what 
happens close to the black hole singularity where quantum gravity effects are dominant and the 
classical Schwarzschild solution is not correct. 
The main results are : 

1. the curvature invariant and the inverse volume operator have a finite spectrum in all the region 
inside the horizon and we can conclude that the classical singularity disappears at the kinematical 
level; on the other side for large eigenvalues of the volume operator we find the classical inverse 
volume behavior, 

2. the solution of the Hamiltonian constraint gives a difference equation for the coefficients of the 
physical states defined in the dual space of some dense subspace of the kinematical Hilbert space. 
All the coefficients in the difference equation are regular in the classical singular point then we 
have a solution of the singularity problem also at the dynamical level. 

An important consequence of the quantization is that, unlike the classical evolution, the quantum 
evolution does not stop at the classical singularity and the "other side" of the singularity corresponds to 
a new domain where the triad reverses its orientation. From the difference equation we obtain physical 
states as combinations of a countable number of components of the form Vw+n/io,^ i+tomo lA'-E + 
n/iQj^E 1 + m Mo) (where ~ 1 at the Plank scale and n,m € Z); any component corresponds 
to a particular value of the volume of the space section. We can interpret He as the time and the 
anisotropy he 1 as the space partial observable [23j that defines the quantum fluctuations around the 
Schwarzschild solution. We recall that \E\ = b 2 , therefore in the classical theory and in quantum 
mechanics, we can regard \E\ as an internal time. So the function VViE+^Mo.p i+m^t i s the wave 
function of the Black Hole inside the horizon at the time he + tiho and we have a natural and 
regular evolution beyond the classical singularity point which is in he = localized. A solution of the 
Hamiltonian constraint corresponds to a linear combination of black hole states for particular values 
of the anisotropy he 1 at the time he- 

It is interesting to recall that beyond the classical singularity the eigenvalue he is negative and so 
we can suggest a new universe was born from the black hole formation process. In LQBH scenario 
pure states which fall into black hole emerge in a new universe as pure states and the information 
loss problem is avoided. Information is not lost in the black hole but it exists again in the space-time 
region in the future of the avoided singularity. 
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